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Abstract 

The work is devoted to the relaxation limit in larger Besov spaces for compress- 
ible Euler equations, which contains previous results in Sobolev spaces and Besov 
spaces with critical regularity. Such an extension depends on a revision of com- 
mutator estimates and an elementary fact which indicates the connection between 
homogeneous and inhomogeneous Chemin-Lerner spaces. 
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1 Introduction 

In this paper, we consider the following nondimensional compressible Euler equations 

d t p + V-(pv) = 0, 

dt(pv) + V • (pv ® v) + Vp(p) = -pv/r { ' ' 

for (t,x) £ [0,+oo) x 1^ with N > 1. Here p = p(t, x) is the fluid density function; v = 
(v 1 , v 2 , v d ) T (T represents the transpose) denotes the fluid velocity. The pressure p{p) satisfies 
the classical assumption 

p'(p) > 0, for any p > 0. 
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An usual simplicity p(p) := /9 7 (7 > 1), where the adiabatic exponent 7 > 1 corresponds to 
the isentropic flow and 7 = 1 corresponds to the isothermal flow. The nondimensional number 
< r < 1 is a (small) relaxation time. The notation V, <8> are the gradient operator (in x) and 
the symbol for the tensor products of two vectors, respectively. 
System (jl.ip is complemented by the initial conditions 

(p,v)(0,x) = (po,v ). (1.2) 

For fixed r > 0, as we all know, the relaxation term which plays the role of damping can 
prevent the finite time blow-up and the Cauchy problem (|1.1|) - (|1.2|) admits a unique global 
classical solution, provided the initial data is small under certain norms. In this direction, 
such problem was widely studied by many authors, see e.g. (HJ EJ EH 151 EU E2] [25] ES] and 
references therein. In addition, it is proved that the solutions in [25] has the L°° convergence 
rate (1 + i) _3//2 (iV = 3) to the constant background state and the optimal L p (l < p < 00) 
convergence rate (1 + i) _ ^/ 2 ( 1_1 /p) in general several dimensions [26j . respectively. In one space 
dimension in Lagrangian coordinates, Nishida [21] obtained the global classical solutions with 
small data, and the solutions following Darcy's law asymptotically as time tends to infinity was 
shown by Hsiao and Liu [8]. For the large-time behavior of solutions with vacuum, the reader 
is referred to [9j [10] . Nishihara and Yang [22] studied the boundary effect on the asymptotic 
behavior of the solutions to the one-dimensional initial-boundary value problem. Later, Liu and 
Wang [15] considered the 2-D initial-boundary value problem in the wedge-space, and proved 
the global existence of classical solutions. 

Another interesting line of research is to justify the singular limit as r — > in (II. ip . To do 
this, we introduce the time variable by considering an "(D(1/t)" time scale 

(p r ,v r ){ 8 ,x) = (p,v)(^). (1.3) 

(1.4) 



Then 



^ + V-(^) = 0, 

7*a.(£f-) + r 2 v • (^T^) + £f- = - w 



with the initial data 



(p T ,v T )(x,0) = ( /3o ,vo). (1.5) 



At the formal level, at least, if we can assume that p — is uniformly bounded, it will be shown 
that the limit N of p T as r — > satisfies the classical porous medium equation 

d s M - Ap(M) = 0, 

Af(x,0) = p , li ' Dj 

which is a parabolic equation since p{M) is strictly increasing. 

This singular limit from hyperbolic relaxation to parabolic equations have attracted much 
attention, see [5j EJ HHJ H7J HHJ [20], [27] and therein references. The relaxation results mentioned 
for smooth solutions fell in the framework of the classical existence theory of Kato and Majda[TT] 
116] . The regularity index of Sobolev spaces (in x) H a (M. ) is assumed to be high (cr > 1 + N/2 
with integer) . Recently, the first author and Wang [29] studied the limit case of regularity index 
(<7 = 1 + N/2) where the classical theory fails. They developed a new commutator estimate to 
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overcome the technique difficulty and constructed global classical solutions in the critical Besov 
spaces f?2~i 2 (R JV ). Furthermore, based on the Aubin-Lions compactness lemma in [23], it was 
justified that the (scaled) density converges to the solution of the porous medium equation. 

The main purpose of this paper is to generalize the relaxation limits in Sobolev spaces with 
higher regularity and Besov spaces with critical regularity. Due to the elementary fact that 
Sobolev spaces H a (K ) := B^ili ), a natural question is whether those results hold in larger 
Besov spaces B% r {R N ) or not, whose indices satisfy the following condition: 

a > 1 + N/2, 1 < r < 2 or cr = l + N/2, r = 1. (1.7) 

In the present paper, we shall answer the question. The main difficulty lies in the a priori 
nonlinear estimates, in particular, the commutator estimates. To overcome it, we need a more 
general version of commutator estimates in Proposition 12.31 (also see [7]), which relaxes the 
restriction on the couple (s,r). For completeness and the reader convenience, we present the 
proof in Appendix with the aid of the Bony's decomposition. Besides, this extension also heavily 
depends on an elementary fact developed for general hyperbolic system of balance laws (see 
Lemma 14.11 or [28]), which indicates the connection between homogeneous and inhomogeneous 
Chemin-Lerner spaces. Precisely, our results are stated as follows. 

Theorem 1.1. Let the couple {o~,r) satisfy the condition \1. 7| ) and let p > be a constant 
reference density. There exists a positive constant 5q independent of r such that if 

||(po-p,m )|| B ^ (H ,jv) < 5 

with mo := po v o> then the Cauchy problem U. 1\) - IT3$I has a unique global classical solution 
(p,m) € C 1 (R+ x R N ) satisfying (p — p, m) G C(B^ r (R N )) n C 1 (B^~ 1 (R N )). Furthermore, the 
global solutions satisfy the inequality 



||(p-p,m)|| £oo(B? rRjv ^ +^ ( — ro/ _ ,„ w „ + V^(Vp,Vm) 



2.r 



L 2 (B. 



CT-l 

2.r 



< CoIKpo -P,mo)llBS ir (R"), ( L8 ) 
where m := pv, po and Cq are some uniform positive constants independent o/r(0 < r < 1). 



Remark 1.1. It should be pointed out Theorem 11.11 contains previous results (see [T4"l 
[26| [29] and the references therein) for compressible Euler equations. The proof depends on the 
Fourier localization methods and Shizuta-Kawashima algebraic condition which has been well 
developed by the second author et al. [I2j [30] for generally hyperbolic systems of balance laws. 
In particular, a concrete entropy for current compressible Euler equations is available, which 
has been verified in, e.g., |14[ 128]. In addition, we track the singular parameter r in the uniform 
energy inequality (jl.8p . which plays a key role in the analysis of relaxation limit problem. 

As a direct consequence, using the standard weak convergence method and Aubin-Lions 
compactness lemma (|23J), we further have the analogue relaxation convergence as in |29] , 

Theorem 1.2. Let (p, m) be the global solution of Theorem \l.l[ Then it yields 
p T — p is uniformly bounded in C(1E> + , B% r (W N )); 
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is uniformly bounded in L 2 (R + , B% r (R N )). 

T ' 

Furthermore, there exists some function J\f G C(R + , fi + B% r {M. N )) which is a unique solution of 
&1.6\) . For any < T, R < 00, {p T (s,x)} strongly converges to AT(s,x) in C([0,T], (B%~ S (B r )) 
as t — > 0, where 5 E (0, 1) and B r denotes the ball of radius r in R . In addition, it holds that 

||(7V(s, -) -p\\ B « r{RN) < Co\\{po - p,m )\\ Blr{RN) , s > 0, (1.9) 

where Cq > is a uniform constant independent of r . 

Remark 1.2. Theorem 11.21 gives a rigorous description in larger spaces that the porous medium 
equation is usually regarded as an appropriate model for compressible inviscid fluids in small 
amplitude regime of relaxation time r. In comparison with our recent results in [29], Theorem 
I1.HI1.2I also hold in the cases of general pressure and arbitrary space dimensions except for the 
regularity consideration. 

Remark 1.3. From the symmetrization in Sectj3j we see that the dependence of the matrices 
A J (j = 0, 1,2 • --,N) with respect to the total variable W rather than only W2. The concrete 
context of matrices enables us to obtain the uniform frequency-localization estimates and the 
relaxation limit. However, to the best of our knowledge, it is unknown for generally hyperbolic 
systems to get corresponding results. Therefore, this paper can be regarded as the effort to the 
open question in [3] (Remark 15, p. 225). 

The rest of this paper unfolds as follows. In Sect. El we briefly review the Littlewood- 
Paley decomposition and properties of Besov spaces and Chemin-Lerner spaces. In Sect. El we 
reformulate the equations (jl.ip as a symmetric hyperbolic form in terms of entropy variables, 
and the local existence and blow-up criterion of classical solutions in critical spaces are presented. 
Sect. H] is devoted to deduce the a priori estimates in Chemin-Lerner spaces by using Fourier- 
localization arguments, which are used to achieve the global existence of classical solutions. 

2 Preliminary 

Throughout the paper, / < g denotes / < Cg, where C > is a generic constant, f ~ g 
means / < g and g < /. Denote by C([0,T], X) (resp., C 1 ([0,T],X)) the space of continuous 
(resp., continuously differentiable) functions on [0,T] with values in a Banach space X. Also, 
||(/, g, h)\\x means H/Hx + where /, g G X. (f,g) denotes the inner product of two 

functions f,g in L 2 (R N ). 

In this section, we briefly review the Littlewood-Paley decomposition and some properties 
of Besov spaces. The reader is also referred to, e.g., [1\ for more details. We omit the space 
dependence, since all functional spaces (in x) are considered in M. N . 

Let us start with the Fourier transform. The Fourier transform / of a L 1 -function / is given 

by 

Ff= j f(x)e- 2 ™<dx. 

More generally, the Fourier transform of any / G S', the space of tempered distributions, is 
given by 

{Ff,g) = (f,Fg) 
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for any g G S, the Schwartz class. 
First, we fix some notation. 

S = [<pe S,d a Ff{0) = 0,Va G N w multi-index}. 

Its dual is given by 

S' = S'/V, 

where V is the space of polynomials. 

We now introduce a dyadic partition of M. N . We choose 4>o G S such that </>o is even, 

supp^o := A = |e G R N : - < \£\ < jj}, and 4> > on A . 
Set A g = 2 q Ao for g£Z. Furthermore, we define 

and define f g 6 5 by 

».(0 - 14,(0 



T, qe z MO' 

It follows that both J 7 ^^) and $ ? are even and satisfy the following properties: 
and 

As a consequence, for any / G Sq, we have 

oo 

E ***/ = /■ 

g=— oo 

To define the homogeneous Besov spaces, we set 

A J = $,*/, g = 0,±l,±2,... 
Definition 2.1. For s£K anc? 1 < p, r < oo, £/te homogeneous Besov spaces Bp r is defined by 

B;, r = {f£S> : \\f\\ Hr < oo}, 

where 



B 3 ~ ' 



E 9e z(2 9S ||Aj||L P r) 1/r , r<oo, 
, sup (?gZ 2'i ,s ||A (? /|| L P, r = oo. 
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To define the inhomogeneous Besov spaces, we set ^ G Cq°(R ) be even and satisfy 

oo 

^*(o = i-E^(o- 

It is clear that for any / € S', yields 

oo 

We further set 

f 0, g < -2, 

A,/=^**/, g = -l, 

($ q *f, q = 0,1,2,... 

Definition 2.2. For s£l and 1 < p,r < oo, the inhomogeneous Besov spaces B pr is defined 
by 

B' p , = {f€S': \\f\\ Bkr < oo}, 

where 

Bs = /(E^-i(2"||A,/|W) 1/r > r<oo, 
lsup,>_ 1 2«»||A,/|| L ,., r = oo. 

Next we turn to Bernstein inequalities. 

Lemma 2.1. Lei G N and < i?i < i?2- There exists a constant C, depending only on R\,R2 
and N , such that for all 1 < a < 6 < oo and / G L a , 

SuppF/ c{(eR N : |e| < R\\} sup ||a a /|| L6 < C k+1 X k+d ^\\f\\ La ; 

\a\=k 

SupjxFf C^el^: i?iA < 1^1 < i? 2 A} C'^X^lfUa < sup ||d a /|| L a < C fc+1 A fc ||/|| La . 

\a\=k 

As a direct corollary of the above inequalities, it holds that 
Remark 2.1. 
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C B pr 



+ M < IF/II^ <C||/|L S+H ; 



l^/lli?^ < C||j || B »+|a|, 



for all multi- index a. 

The Besov spaces defined above obey various inclusion relations. In particular, we have 
Lemma 2.2. Let s£R and 1 < p, r < oo, then 

(1) lfs>0, then B s pr = L p C\ B s p y, 

(2) Ifs<s, then B s p;r ^ B; -; 
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(3) B»{ p ^ Co, B%* ^ C (l < V < oo); 

where Co is t/ie space 0/ continuous bounded functions which decay at infinity. 

On the other hand, we also present the definition of Chemin-Lerner space-time spaces ini- 
tialed by J.-Y. Chemin and N. Lerner [2], which are the refinement of the spaces Lj,(Bp r ) or 
L T (Bp r ). 

Definition 2.3. For T > 0, s G R, 1 < r, < 00, f/ie homogeneous mixed time-space Besov 
spaces L e T (Bp r ) is defined by 

l 6 T (B^ r ) := {/ G L>,T;S ) = ll/llzg,(B' r ) < 

where 

:=(E( 29 l A ^H^)) r )' 

with the usual convention if r = 00. 

Definition 2.4. For T > 0,s £ < r,9 < 00, the inhomogeneous mixed time-space Besov 
spaces L^(Bp r ) is defined by 

L e T (B s p , r ) ■= {/ G L>,T;S') : ll/llz^,.) < +~}, 

11/11^ p) :=(E( 29S iiN/ii^(L P) r)" 
?>-i 

with the usual convention if r = 00. 

We only state some basic properties on the inhomogeneous Chemin-Lerner spaces, since the 
similar ones follow in the homogeneous Chemin-Lerner spaces. 

The first one is that Lfp{B pr ) may be linked with the classical spaces L^,(Bp r ) via the 
Minkowski's inequality: 

Remark 2.2. It holds that 

Let us also recall the property of continuity for product in Chemin-Lerner spaces L^p{B s pr ). 
Proposition 2.1. The following inequality holds: 

ii^iiz§,(fl- P ) ^ c ^fk e T HL^)\\9\\zp {B s r) + yw^^yh^^sj 

whenever s > 0, 1 < p < 00, 1 < 9, 9\, 62, #3, #4 < 00 and 

1 _ 1 11 1 
9 9\ 92 #3 #4 

As a direct corollary, one has 

\\fg\\ie T{B s :r) ^cil/ll^i^jNIz^^.,.) 

whenever s > d/p, 3 = g~ + 0~- 
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In the next symmetrization, we meet with some composition functions. The following con- 
tinuity result for compositions is used to estimate them. 

Proposition 2.2. Let s > 0, 1 < p,r,0 < oo, F' G W^ +1 '°°(I; R) with F(0) = 0, T G (0,oo] 
and u G L 6 T {B s p r ) n L~(L°°). T/ien 

IW)IIzs.(b* p ) ^ ^ + Il/ll^(^)) w+1 ll^ll^.-ll/llz5. (fl . ir) - 

In addition, we present some estimates of commutators in homogeneous and inhomogeneous 
Chemin-Lerner spaces to bound commutators. 

Proposition 2.3. /?/ Let 1 < p < oo, 1 < < oo and s G (— ^ — 1, ^]. T/ien i/jere exists a 
generic constant C > depending only on s,N such that 

( ||[/,A 9 ]< 7 || LP <Cc (? 2-^+ 1 )||/||. i v +1 ||< 7 ||^ 

' ||[/ ! A ? ] 5 || i , ( ^ ) <c C9 2-^+ 1 )||/|| v +1 \\g\\ z e 2{ ^ , 
' 4 ( B P i ) T p,x 

with 1/0 = l/6\ + 1/02, where the commutator [•,•] is defined by [f,g] = fg — gf and {c q } 
denotes a sequence such that ||(cg)||ji < 1. 



3 Entropy and Local well-posedness 

First, let us introduce the energy function which is just an entropy in the sense of Definition 2.1 
of TI\ : 



Iml 2 f p p'is^ 

ri(p,m) := h hip) with m = pv and h'ip) = 

2p J 1 



ds. 



s 

For this rigorous verification, see, e.g., [28]. Furthermore, the associated entropy flux is 

'ml 2 , , . \ m 



Define 



Clearly, the mapping U — > W is a diffeomorphism from C( Pim ) := M + x M. d onto its range 0\y, 
and for classical solutions (p, v) away from vacuum, is equivalent to the symmetric system 

d 

A°iW)W t + Y^A'(W)W Xj = H(W) (3.1) 
j'=i 

with 

A°iW)=( 1 ^ 

1 ; \W 2 W 2 ®W 2 +p'{p)I d 



AHW)=[ W * W?W 2j +p>(p)e] 

' 1 ' 1 W 2 W 2j +p'{p)e j W 2j {W 2 ® W 2 +p'(p)I d )+p'{p)(W 2 ®e j + e j ® W 2 ) 



H(W) = G(U(W))=p'(p)( _^ 



where Id stands for the dxd unit matrix, and ej is (i-dimensional vector where the jth component 
is one, others are zero. It follows from the definition of entropy variable W that h'(p) = 
W\ + | W 2 \ 2 /2, so p'{p) can be viewed as a function of W, since p is the function of W± + \W 2 \ 2 /2, 
i.e. ofW. 

The corresponding initial data become into 

I 1 2 

W(0,x) :=W = (-^+/ i '(p ),vo). (3.2) 

In [28], we have recently established a local well-posedness theory for generally symmetric 
hyperbolic systems in the framework of critical spaces, which is regarded as the generalization 
of the classical existence theory of Kato and Majda 116] . Actually, the theory is also true 
in larger Besov spaces, whose the regularity indices satisfy the condition (jl.7p . Hence, we can 
have the following local existence result for the concrete problem (I3,lj) - (|3.2I) , 

Proposition 3.1. For any fixed relaxation time r > 0, assume that the initial data W$ satisfy 
Wo — WE. B^jXW := (h'(p),0)) and take values in a compact subset ofOw, then, there exists 
a time To > such that 

(i) Existence: the Cauchy problem \3.1\) - (3lfy has a unique classical solution W G C 1 ([0,To] x 
M. d ) satisfying 

W-WeOr^B^n^B^ 1 ); 



(ii) Blow-up criterion: there exists a constant Co > such that the maximal time T* of 
existence of such 

T* is finite, then 



existence of such a solution can be bounded from below by T* > -7^? — Sr< • Moreover, if 

j y - \\w -w\\ B <y ' J 

2,r 



limsup \\W — W\\b% = oo 

if and only if 

r T* 

\\VW\\L°°dt = oo. 



4 Global well-posedness 

This Section is devoted to the global existence result in Theorem ll.il To show that the solutions 
of (|3.ip - (|3.2p . are globally defined, we need further a priori estimates. 

To do this, for any time T > and for any solution W — W G Ct{B 2 ^ H C^B^' 1 ), we 
define by E(T) the energy functional and by D T (T) the corresponding dissipation functional: 

E(T) := \\W-W\\ Z ( j, 
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D T (T) -^ll^llz^^ + ^llV^llza^), 

and E(0) := \\W - W\\ B - r - We also define 

S(T) := ^^^([o^jxRd) + ||VW|| iO o([ 0]T ] xR d). 

Note that the imbedding in Lemma 12.21 and Remark I2.2( we have S(T) < CE{T) for some 
generic constant C > 0. 

The next central task is to construct the desired a priori estimate, which is included in the 
following proposition. 

Proposition 4.1. Let W be the solution of (E2P-(ED satisfying W-W G C T {B^ r )nC^(B^ x ). 
IfW(t,x) takes values in a neighborhood of W , then there exists a non- decreasing continuous 
function C : IR + — > M + which is independent of t, such that the following nonlinear inequality 
holds: 

E(T) + D T {T) < C(S(T))(e(0) + E{T) 1/2 D T {T) + E(T)D T (T)^ . (4.1) 

Furthermore, there exist some positive constants 5i,(ii and C\ independent of t, if E{T) < 6%, 
then 

S(T)+Mii?r(r)<CiS(0). (4.2) 

The proof of Proposition 14.11 in fact, is to capture the dissipation rates of W = (Wi,W2) 
in turn by using the Fourier localization arguments. First, we give an elementary fact that 
indicates the connection between the homogeneous and inhomogeneous Chemin-Lerner spaces, 
which have been established in the recent work [28]. Precisely, it reads as follows: 

Lemma 4.1. Let s > 0, 1 < 9,p, r < +oo. When 9 >r, it holds that 

LULnnl e T (B^ r ) = L 9 T (B^) (4.3) 

for any T > 0. 

Next we begin to prove the Proposition 14.11 For clarity, we divide it into several lemmas. 

Lemma 4.2. Under the assumptions stated in Proposition \4-l[ there exists a non- decreasing 
continuous function C : M + — > M + which is independent of t, such that the following estimate 
holds: 

\\ W -W\\l~(bz ) + ^W W 4lUb ? ) <C{S{T)){E{Q)+E{T) l / 2 D T {T)). (4.4) 
Proof. The proof is divided into three steps. 

Step 1. The Lf{B^ r ) estimate ofW and the L\{B^ r ) one of W 2 
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Applying the homogeneous localization operator A q (q E Z) to (|3.ip . and taking the inner 
product with A q , we obtain 

d 

A°(W)A q W t + Y J A j {W)A q W X] + LA q W 

i=i 

d 

= -[A q ,A°(W)]d t W -J2l&<i,^(W)}d Xj W + r(W) (4.5) 

i=i 



with 

L ^ W '= ~r ( P>(Aw 2 ) aM ^ : = \ ( A 9 [(^(P)-P'(P))^] 



where the commutator [-, •] is defined by [/, g] := fg — gf. 
Take the L 2 -inner product of (j4.5|) with A q W to get 



1 (A\W)A q W,A q W)t + ^\\A q W 2 f L 2=y / Ii(t), (4.6) 



t=l 

where 



/i := i(^°(W)A,W; A g W), J 2 := ~ J2{d Xj A>(W)A q W, A q W), 

3=1 

d 

J 3 := -([A,,A (W)]W,AgWO, I 4 ^-^([Ag.A^W^W.A^), J s := (r(W),A ff W>. 

3=1 

In what follows, we begin to bound these nonlinear energy terms. Firstly, remark that 
A°(W) = ( 1 W ? ) + (° 

y > V w 2 w 2 ® w 2 J V pW* 

:= A I (W 2 ) + A° II (W), (4.7) 



we further get 



n i"(ir i - ( 9 ' W 2 T ) ( i - 

a * nKK2; 1 d t w 2 d t w 2 ®w 2 + w 2 ®d t w 2 I ' " 



and 



= ( o swDrZmwu ) • < 4 - 9 ' 

where Dyy stands for the (row) gradient operator with respect to W . From the decomposition 
in (14, 7ft of A (W), I\ can be written as the sum 

h(t) = ^(d t A j(W 2 )A q W,A q W) + ^(d t A n (W)A q W,A q W) =: I u (t) + I l2 (t). (4.10) 
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With the aid of the embedding properties in Lemma 12.21 and Remark 12.11 we obtain 



T 



hi(t)dt 



(d t w 2 ■ A,w 2 A 9 wi + -\A q w 2 \ 2 w 2 ■ d t w 2 ] 

l A 



< C(5(T)) / \\d t W 2 \\ Loo \\A q W 2 \\ L 2\\A q W\\ L 2dt 



VW\\ L oo + -||W 2 ||ioo \\A q W 2 \\ L 2\\A q W\\ L 2dt 



< C(S(T)) 

+ -||A ? W|| i o ?(i 2 ) ||W 2 || i 2 (Loo) ||A 9 W 2 || L 2, (i 2 ) . 

Furthermore, multiplying the factor 2 2qa on both sides of (14. IIP gives 



(4.11) 



r 



I xl {t)dt < C(5(T))^||W-W|| £ . ( ^ r) (||^ 2 || z , (B2%) ||VW|| z , ( ^ ;1) 



+-|I^|| £|( ^ ;1) II^ 2 |Iz|(b ? j 



(4.12) 



Here and below, {c q } denotes some sequence which satisfies ||(c g )||/r < 1 although each {c q } is 
possibly different in (|4.12p . Allow us to abuse the notation for simplicity. Similarly, we have 



[ T ' Iw(t)dt < C(S(T))2 2qr7 f T [ \d t W\\A q W 2 \ 2 dxdt 

JO JO JR d 



< C(S(T))2 2 ^( \\VW\\ L ~ iLoo) + -\\W 2 \\ L 



J (L°°) 



\ a q w ^\\l 2 t (l^ 



< c(s(r))4i|w-w||^, ( ^ p) -||w a ||^ (V 



(4.13) 



where we used the smallness of r(0 < r < 1). 

Combining (|4.12p - (|4.13p . it follows from the basic fact in Lemma 14. II and Young's inequality 
that 



( T h(t)dt <C(S(T))c 2 q E(T)D T (T) 2 . 
Jo 



(4.14) 



Secondly, we turn to estimate I 2 . For this purpose, we write 



A j (W) 



W 2j W?W 2j 



+ 



P'(p)e] 



W 2 W 2j W 2j {W 2 ®W 2 ) 


p'(p)ej W 2jP '{p)I d + p'{p) (W 2 ej + ej W 2 ) 
:= A{{W 2 ) + A J n {W). 

Therefore, it is not difficult to get 

cT i d 

2 



(4.15) 



1 d 

- J2(d Xj A 3 I (W 2 )A q W, A q W)dt 



3=1 
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< C(S(r))2 2 ^||VW 2 || L|(Loo) ||A 9 ^|| L? c (L2) ||A^|| L 2 (i2) 

< C(S(T))c 2 q \\W - mz nBl jW2\\~ LHB , r) 



< C(S(T))4\\W-W\\~ L¥{Blr)[ - 

< C{S{T))c 2 q E{T)D T (T) 2 . 



-ii^iiW) +r||w|l W) 



(4.16) 



Thanks to the null block located at the first line and the first column of A^iW), we are led to 
the estimate 



1 d 

- Y / (9 Xj A 1 II (W 2 )A q W, A q W)dt 



3=1 



< C(S(T))2 2qa [ \\VW\\ L °o\\A q W 2 \\ L 2\\A q W\\ L 2 
Jo 



dt 



< cr(5(r))^||w-w||^ w (-||w 2 ||^ w 

< C{S{T))c 2 q E{T)D T (T) 2 . 
Hence, together with (I4.16| )-( l4.17jl . we arrive at 



^qcr 



hdt 



< C(S(T))c 2 E{T)D T (T) 2 . 



(4.17) 



(4.18) 



Thirdly, we also use the decomposition (|4.7p of ^4°(VF) in order to estimate the commutator 
occurring in 1%. Using the commutator estimate f)5.2f) in Proposition 15-H we get 



■([A q ,A%W 2 )}d t W,A q W)dt 



< 2 2 ^\\[A q ,A%w 2 )]d t w\\ LlT{L2) \\A q w\\ L¥(L2) 

+IIWIl24(L«)||A?(w 2 )||^ ( ^J||A,w|Uc ?(£3) 

< C2«°c q \\A°{W 2 )\\ l2T{ ^ 



< C(S(T))c 2 q \\W-W\\ l¥{K) 

x (llv^ll z , (B2V) + i|r 2 || 4(B2V) 

< C(S{T))c 2 q E{T)D T (T) 2 , 



(4.19) 



where (a, r) satisfies the condition (jl.7p and we take 6\ = 6 2 = 63 = 6^ = 2 in (|5.2p . 

Recall that p'(p) can be viewed as a function of W and the nice construction of A® n (W), we 
gather that 



c^lqa 



■([A q ,A° n (W)]d t W,A q W)dt 
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< 2 2qa 



[k q , P '{ P )I d }d t W 2 \\ L 2\\K q W 2 \\ L 2dt 



< C2^c q (\\W(p)\\L ¥{L ^)\\dtW 2 \\^ rrl) + \\d t W 2 \\ L 2 iLoa) \\p'(p)\\~ L¥{ ^ r , 



< C(S(T))c 2 q E(T)D T {T) 2 , (4.20) 

where we used the commutator estimates in Proposition 15,11 and the homogeneous case of Propo- 
sition [2T21 

Therefore, we finally get 



^qcr 



hdt 



< C(S(T))c 2 E(T)D T (T)' 



(4.21) 



The estimates of commutators arising in I4 are actually simpler than (j4.19|) - (j4.20p . since 
they only involve spatial derivatives. Let us give the simplified steps only: 



c^lqc 



^2{[A q ,A j j(w 2 )]d Xj w,A q w) 



i=i 



dt 



< C(S(T))c 2 E(T)D T (T) 2 , 



(4.22) 



and 



< C2 2qa 



J2([A q ,A j n (w)]d Xj w,A q w) 

3=1 



dt 



T 



< C{S{T))c 2 q \\W-W\\ 7 ^ (m J^2|lz^ ( ^ r) ||VW|| 4( ^ ;1) 



< C(S(T))c 2 E(T)D T (T) 2 , 

where A J II1 (W) and A J II2 (W) can be viewed as the smooth functions of W. 
Hence, from (j4T22l - (|4T23D . we deduce that 



c^lqa 



hdt 



< C(S(T))c 2 E{T)D T (T) 2 . 



Finally, due to the homogeneous cases of Propositions I2.1H2.21 ^5 can be estimated as 



c^lqa 





2 2qa r 


[ hdt 




Jo 





\A q [(p'(-p)- P '{p))w 2 }\\ L 2\\A q w 2 \\ L 2dt 



5. 

T 



< ^il^(^-y(p))W2|| Z|( ^ r) ||W 2 || Z|(B . r) 

1 2 



< C-(5(T))^[|W-^ ( ^ p) [|W 2 [|^ r) 

< C(S(T))c 2 q E(T)D T (T) 2 . 



WlA^A^iWyd^WWAqWWv + \\[A q ,A 3 II2 (W))d x .W\\\\A q W 2 \\ L 2)dt 

'/H " " Hi|?(BJ r )ll " -11 ^2 / R CT 



(4.23) 



(4.24) 



(4.25) 
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In addition, we define 

\ LK := (A°(W)f, f) 1 / 2 . 



'AO 



Clearly, the norm ||/||j^a ~ ||/||l2, since W(t,x) takes values in a neighborhood of W. Conse- 

A° 

quently, according to inequalities (|4.14[) . f|4.18[) . (|4.2ip and (|4.24p - (|4,25|) . we conclude that 

2 2 «°\\A q (W(t) - W)\\\ 2 + — [|A,W 2 ||| ?(ra) 

< C(5(r))2 29CT ||A (J (VFo - W)||| 2 + C(S(T))c 2 £(T)D T (T) 2 . (4.26) 
Then it follows from the classical Young's inequality that 

2^\\K q {w - woiUoo (£a) + ^||A 9 w 2 |L2 (i2) 

< c(5(r))2^||A 9 (Wo-Tr)|| i2 + c(5(r))c g £;(r) 1 /2£) T ( r ). ( 4 . 27 ) 

Hence, taking the £ r -norm on q £ Z on both sides of (|4.27p gives immediately 

11^ " W h ¥ (B lr ) + ^ll^llz^^) < C(S(T))(E(0) + E{Tfl 2 D T (T)). (4.28) 



Step £ TTie L^{L 2 ) estimate of W and the h\(L 2 ) one ofW 2 
It is convenient to introduce the relative entropy fj(p,m) by 

fj(p, m) := r](p, m) - rj(p, 0) - ^ p r/(p, 0)(p - p). (4.29) 

Then the strictly convex quantity fj(p, m) satisfies 

fj(p, m) > 0, fj(p, 0) = 0, D p fj(p, 0) = 0. (4.30) 

Furthermore, fj(U) is equivalent to the quadratic function \p — p| 2 + |m| 2 and hence to \W — W\ 2 , 

since \p — p\ + |m| lies in a bounded set. Accordingly, the entropy flux qj(p, = 1,2. • ••, d) 
are modified as follows 

qj(p,m) = q j (p,m) - D p r)(p,0)mj, (4.31) 
and we get the entropy-entropy flux equation 

d ■ 1 

d t r}(p,m) +^2d Xj q j (p,m) = p|v| 2 . (4.32) 

i=i 

Integrating (]4.32p over [0, T] x M. N implies immediately 

\\W-W\\ L¥{L 2 ) + ^\\W 2 \\ L 2 t{l2) < C(S(T))\\W -W\\ L 2 

< C{S{T))(E(0) + E(T) l ' 2 D T (T)). (4.33) 

Step 3. Combining the above analysis. 

Note that the fact $0]), flHD is followed by IjOBl) and (P3l) . □ 
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The next step consists in deriving the Lj<(B%~ ) estimate of VW. To this end, we use an 
important stability condition ( "Shizuta-Kawashima" condition), which has been first formulated 
in [23] by the second author, and well developed in [12|, I28 |, [30] . As shown by the recent work [28], 
the compressible Euler equations with relaxation satisfies the stability condition at the constant 
state naturally. Furthermore, we obtain the concrete content of the compensating matrix K(£), 
which has been pointed out by [SKH]. Precisely, 

Lemma 4.3 (Shizuta-Kawashima). For all £ G ~R N , £ ^ 0, there exists a (N + 1) x (N + 1) 

matrix if (£) depending smooth on the unit sphere S^ -1 : 



K(0 



1 £ 

1 



such that K(£ t )A°(W) is skew- symmetric and 



f lei o 



K(0j2(jA j (W)=[ 'J' _ pl{p) m ), (4.34) 



3=1 



where A and A 3 are the matrices appearing in the symmetric system \3.1 



Next we shall use the skew-symmetry of the compensating matrix if(£) in the Fourier spaces 
to establish the estimate of VW. 



Lemma 4.4. Under the assumptions stated in Proposition \4-l\ there exists a non- decreasing 
continuous function C : M + —> R + which is independent of t, such that the following estimate 
holds: 

V^ll WlU^-i^ < C(S(T)) (e(0) + E(T) l ' 2 D T {T) + E(T)D T (T)) . (4.35) 



Proof. Firstly, we linearize (|3.ip around the constant state W: 

d 



A°(W)d t W + ^2A j {W)d Xj W = -LW + Q (4.36) 

3=1 



with W :=W — W, where 

t v p (p)w 2 



and 



d 



^2 A°(W)[A°(W)~ 1 A j (W) - A°(W)~ 1 A j (W)}d Xj W 



3=1 



-^A\W)\p'{p)A\W)- 1 -p'(p)A\W)- 1 ] ( ^ Y (4-37) 
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Applying the inhomogeneous localization operator AJq > —1) to (|4.36p implies 

d 

A°(W)d t A q W + Y, Aj (W)d Xj A q W = —LA q W + A q g. (4.38) 
i=i 

Take the Fourier transform (in the space variable x), multiply by — ir(A q W)* K (£)(* transposed 
conjugate), and compute the real part of each term in the resulting equality to get 

rIm((^)*(ir(OA^Tr))^A^+r(A^)*if(e)(^e^W)A^ 

3=1 

= -p'(p)Im((A q Wi)*^Aalv 2 ^ +Tlm((A q W)*K(Z)A q ~gy (4.39) 
Using the skew-symmetry of K(£ t )A°(W), we have 

Im((A^)*(^(C)A°(W))^A^) = ^Im((A^)*(K(OA°(W))A^). (4.40) 
Then, with the help of (|4.34p . the left-hand of (|4.39p is bounded from below by 

^Im((^)*(K(e)A (Ty))^)+Tp'(^)|e||^| 2 -Cr|e||A^y 2 | 2 , (4.41) 

where C > is a generic constant independent of r. By Young's inequality, the right-hand of 
(|4.39p is dominated by 

i - — x= f ... n r 

-rp'm\\A q W\ 2 + ^|A 9 iy 2 | 2 + ^|A^| 2 . (4.42) 

Hence, combine the inequalities (|4.4ip - (|4.42p to deduce that 

l -r P '(p)\i\\A^\ 2 < + ^)\A\W 2 \ 2 + ^\A7G\ 2 

-^ t Im[(A^)%K(OA°(W))A q w). (4.43) 

Multiplying (14.43P by |£| and integrating it over [0, T] x M. N , then by Plancherel's theorem, we 
arrive at 

T\\A q VW\\ 2 Ll{L2) < (\\A q Wo\\ 2 L2 + ||A 9 VW ||| 2 ) + (HA^IIIoc^) + ||A,W[||oo (ia) ) 

+i(||A g iy 2 || 2 -, (L2) + \\A 9 VW 2 \\ 2 ^ m ) +r||A^|| 2 , (i2) , (4.44) 

where we used the uniform boundedness of the matrix K(£)A°(W)(£ ^ 0) and the smallness of 
r(0 < r < 1). 
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By multiplying the factor 2 2q ( a ^ and taking the £ r -norm on q > —1 on both sides of (|4.44|) . 
then extracting the square root of the resulting inequality, we are led to 

V^Wvww^-^ < \\Wo - w\\ Blr + \\w - w\\z~ (Blr) 

+^11^11^) + ^11^112,^-1). (4.45) 
Recalling the definition in (|4.34p of G, with the aid of Propositions 12.1112. 21 we obtain 

\\Gk*(B^) < C(S(T))\\W - mz~(B lr ){~mk HBlr) + \^W\\Z H B^))- (4-46) 
Therefore, we get ultimately 

v^llWII^-i) < C(S(T)){\\W - W\\ Blr + \\W- W\\ Zoo{Blr) + ±\\W 2 \\ l2{Blr) 

+ll w -^ll^ w (^ll^b WlP )+^l vwr b(^))} 

< C(S(T)) (E(0) + E(T) l / 2 D T {T) + E(T)D T (T)^j , (4.47) 

which completes the proof of Lemma 14.41 □ 

Proof of Proposition \4-l\ By combining (|4,4p and ()4.35p . we conclude the nonlinear inequality 
(|4.ip . In addition, the inequality (|4.2p follows from ()4.ip and the a priori assumption E(T) < 5± 
readily. Hence the proof of Proposition 14.11 is finished. □ 

Thanks to the standard boot-strap argument, for instance, see [[19], Theorem 7.1] or the 
outline given in [29], we extend the local-in-time solutions in Proposition l3.1l to the global- in-time 
classical solutions of (|3.ip - (|3.2p . Furthermore, we arrive at Theorem ll.il 

The proof of Theorem ll.2l is just the same as that in [29], we thus feel free to skip the details, 
the interested readers is referred to [29] . 

5 Appendix 

In this section, we present a revision of commutator estimates in Proposition 12.31 and relax the 
restriction on the couple (s,r) in comparison with the commutator estimate in Proposition 12.31 
which enables us to construct the a priori estimates in more general Besov spaces. 

Proposition 5.1. For s > —1, 1 < p < oo and 1 < r < oo, there is a constant C > such that 
|| [/, A q ]g\\ LP < Cc q 2~^ (|| v/Hi- IMI^, + \\g\\ Ln ||/||^+i) (5.1) 

and 

\\[f,A q ]g\\LUL P ) < ^ 9(s+1) (llWII 4l(Lco) ll5ll z « 2(%r) + ll5ll 4 3 (LP1) ll/ll z « 4( ^ lr) ), (5-2) 
where 1/p = 1/pi + l/p2> 1/0 = l/$i + l/#2 = l/#3 + l/#4 and c q denotes a sequence such that 

\\{c q )hr<l. 
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Proof. We decompose 

[f,A q ]g:=K 1 + K 2 + K 3 + K i + K 5 (5.3) 

with 

K\ = [Tf, A q ]g, K 2 = R(f,A q g) K 3 = -A q R(f,g), 
K * = T A qg f> K 5 = -A q T g f, 
where T and R stand for the paraproduct and remainder operators (see [2] by J.-M. Bony), which 
are given by T f g = E 9 '<q- 2 A g '/A 3 # = E gG z S q-if\g and R(f, g) = E|g- g '|<i \f\'9- (EH 
is obtained after noticing that 

fg = T f g + T g f + R(f,g). 
From the definition of A q , we have the almost orthogonal properties: 

A p A q f = if \p - q\ > 2, 

A q (S q ^fA p g) = if |p-g|>5. 

Now, we go back to the proof of inequality (|5.ip . For K%, it follows from the Taloy's formula 
of first order and Young's inequality that 

\\KxWlp < C Y, 2 ~ q \\Sk-iVf\\ L ~\\A k g\\ LP 

\q-k\<4 

< cv-^ +1 )||v/|MMI^ • (5-4) 

p,r 

For i^2, we have 

^ = E E (A*/)(Afc,A^) 

k&L |fc-Jf|<l 

= E E (A fc /)(A fe ,A^). 

|fc-<j|<2 |fc-fc'|<l 

Then by Holder inequalities, we obtain 

\\K 2 \\lp < C E E l|Afc/||Lft||A^A^[|LPi 
\k-q\<2 |fc-fc'|<l 

\k-q\<2 

< Cc^-^ll/H^+ilblUpx, (5.5) 

where 1/p = 1/pi + l/p 2 - 
For K 3 , we have 

K 3 = -A ? (E E A fc /A^ 5 ) 

kez\k-k'\<i 

= - E E A,(AA/Afc, 5 ) 

max(k,k')>q-2 \k-k'\<l 
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Using the Holder inequality and Young's inequality for sequences, we proceed K 3 as follows: 
\\K 3 \\ L p < C J2 \\&kf\\LP2\\A k ,g\\ LP1 

max(k,k')>q-2 |fc-fc'|<l 

< C\\g\\ LPl W^kfW 

k>q-3 

< C Cg 2-^ +1 )||/|U s+1 || 5 || iPl , (5.6) 



where s + 1 > is required. 
For we have 



K * = - E S k ^A q gA k f. 

k>q+l 



In a similar manner as K3, we are led to the estimate 



H-FQHlp < C ^ ||A fc /|| L P 2 ||fif||LPi 

k>q+l 

< Cc q 2-^\\f\\ BS+1 \\g\\ LP1 . (5.7) 
For the estimate of K$, we recall that 

K 5 = - J2 & q (Sk-igA k f), 

\q-k\<4 

furthermore, we get 

H-FQIIlp < C Y \\^kf\\LP2\\g\\LPi 
|g-fc|<4 

< C C(? 2-^ +1 )||/||^+i|| 5 || iP1 . (5.8) 

Together with (j5.3[) - ()5.8|) . the inequality (|5.ip follows immediately. 

In addition, the similar process enables us to obtain (|5.2p . whereas the the time exponent 6 
behaves according to the Holder inequality. □ 

Following from the above proof, it is not difficult to obtain the commutator estimates in the 
inhomogeneous case. 

Corollary 5.1. For s > —1, 1 < p < 00 and 1 < r < 00, there is a constant C > such that 
\\[f,A q ]g\\u < Cc q 2-^(\\Vf\\ L ~\\g\\ B s r + \\g\\ LP1 WfW^i) (5.9) 

and 

\\[f,A q } 9 \\LUL*) < ^ 9(s+1) (ll v /ll4 1( ^) ll3ll 4 2 (^) + 11^11^(^)11/11^^))' ( 5 - 10 ) 

where 1/p = l/pi + I/P2, 1/0 = l/0i + 1/02 = 1/03 + 1/04 and c q denotes a sequence such that 
\\{c q )\W<l. 
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